Constructing all entanglement witnesses from density matrices 
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We demonstrate a general procedure to construct entanglement witnesses for any entangled state. 
This procedure is based on the trace inequality and a general form of entanglement witnesses, which 
is in the form W — p — c p I , where p is a density matrix, c p is a non-negative number related to 
p, and I is the identity matrix. The general form of entanglement witnesses is deduced from Choi- 
Jamiolkowski isomorphism, that can be reinterpreted as that all quantum states can be obtained by 
a maximally quantum entangled state pass through certain completely positive maps. Furthermore, 
we provide the necessary and sufficient condition of the entanglement witness W = p — c p I in 
operation, as well as in theory. 
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Quantum entanglement, which is applied to various 
types of quantum information processing such as quan- 
tum computation 0], quantum dense codingQ, quantum 
teleportationQ, quantum cryptography Q, etc., has been 
incorporated as a central notion in quantum informa- 
tion theory It is well known that entanglement can 
be identified by applying all positive but not completely 
positive (PNCP) maps to a given stated, 0. However .it 
is not easy to find and physically realize PNCP maps[8|. 
An equivalent approach of identifying entanglement is 
based on entanglement witnesses (EWs) 0, . EWs are 
observables that completely characterize separable (not 
entangled) states and allow us to detect entanglement 
physically This make EWs one of the main methods 
of physically detecting entanglement. 

Constructing the EW for an entangled state is a dif- 
ficult task, and the determination of EWs for all entan- 
gled states is a nondeterministic polynomial-time (NP) 
hard problem 0- 11 1. Concerning this topic, much work 
has been done for constructing special EWs (for exam- 
ple, Refs. 

0, EMI)- In this Brief Report, we show a 
general form of EWs from density matrices. This general 
form of entanglement witnesses is deduced from a known 
relation; that is, any quantum state can actually be gen- 
erated from a maximally entangled quantum state with a 
completely positive map. The trace inequality indicates 
an EW can be built up as commuting with a given en- 
tangled state. Therefore, we provide a general procedure 
for detecting entangled quantum states. Furthermore, 
we provide the necessary and sufficient condition of the 
entanglement witness W = p—c p I both in operation and 
in theory. 

For our purpose, we first consider quantum states on 
the finite dimensional Hilbert space Hab = Ha <8> Hb- 
Let dim('H J 4) = d-A, dim(Ks) = d-B and dim("H J 4s) = 
(Iab- Denote P+ as the density matrix of the maximally 
quantum entangled state \p) = d A 1 J2i I*) on Ha® 
Ha, where {l*)}^^ 1 are computational bases in Ha, i-e. 



p+ = \m\- 

The relation between any quantum state and a max- 
imally entangled state. Quantum entanglement, which 
is a fascinating feature of quantum theory, underlines 
the intrinsic order of statistical relations between sub- 
systems of a compound quantum system (l4|. In the fol- 
lowing, we show a relation between all quantum states 
and a maximally entangled state. This relation comes 
from a well-known feature called "channel-state duality" 
or "Jamiolkowski isomorphism" or "Choi-Jamiolkowski 



15] 



isomorphism" . This result first appeared in Ref. 
with a proof. We now show it in a different manner. 

Lemma 1. Any matrix H on Hab is a Hermitian 
matrix if and only if it can always be written as 



fr = (z®*)(p+) 



(i) 



where <I> is a hermiticity-preserving linear map. 

Lemma 2. ( t 16]) A linear map A: Ha — > Hb is 
completely positive if and only if the matrix D € Hab 
given by 



D = (/®A)(P+) 



(2) 



is positive semidefinite. 

By Lemmas 1 and 2, we can get the following result. 

Theorem 1. Any matrix p on Hab is a bipartite 
density matrix if and only if it can always be written as 



p=(I®A)(P+), 



(3) 



where A: Ha ~^Hb is a completely positive map. 

Note that A may not be trace preserving and p may 
not be normalized. A similar result on Ha ® Ha was 
shown by DiVincenzo et al. 17]. If A is trace preserving, 
A is an entan glem ent-breaking channel (EBC)[l8j], but 
not vice versa [191 ]. 

The general form of entanglement witnesses. Starting 
from the positive map [a], the concept of EW was ap- 
plied to detecting the presence of entanglement [2pj] . EWs 



2 



are observables whose expectation value can reveal some- 
thing about the entanglement in a given state 21] . A 
Hermitian matrix W = W* on Hab is an EW if it has 
(i) at least one negative eigenvalue and (ii) nonnegative 
mean values in all separable quantum states, or equiva- 
lently satisfy 



in Ua, {mto 1 inHsfll. We have I®pk s (P+) 



By Eqs. 0) and ©, we have 



W = p-p{I® A S )(P+) = p 



c p I. 



pi 



(10) 



(ij,avb\W\havb) > 



(4) 



for all pure product states HI H^- If we have 

negative mean value in a quantum state for an EW, the 
quantum state is entangled. In that case, we say that 
the EW "witnesses" (detects) the quantum state. To 
balance out the "not trace-preserving" property of the 
completely positive map in this Brief Report, we need 
another property of EWs: (hi) if W is an EW, keeps 
all properties of W as an EW for a non-negative number 
7. Note that the third property is different from the 
definition by Lewenstein et al. for comparing the action 
of different EWs0. 

Theorem 2. Any bipartite density matrix tt is entan- 
gled if and only if there exists a density matrix p and a 
non-negative number c p such that the matrix 



W 



c p I 



(5) 



satisfies tr(M / 7r) < and tr(Wa) > for all separable 
states a. 

This result shows that every entangled state in a com- 
posite system has an EW in the simple form W — p—c p I, 
where c p is a non-negative number and p is a density ma- 
trix. This result also shows that the research on density 
matrices can replace the research on EWs since c p I is 
simple. 

Proof. By Lemma 1, any EW W' can be written as 



W' = {I®Q){P+), 



(6) 



where is a positive map [20, |24{. By property (hi) of 
EWs, 



w = {i-p)W = (i-p)(/®e)(p+) 



(7) 



is the same EW as W for < p < 1. 

We could mix (1 — p)0 with a simple completely posi- 
tive map pA s : (1 — p)Q +pA s . By structural completely 
positive approximation (SCPA) and structural physical 
approximation (SPA) [2H 0] for proper p, 

[I ® ((1 - p)0 + pA s )](P+) = I ® A'(P+) = p, (8) 

where A' = (1 — p)0 + pA s is a completely positive map 
and p is a density matrix by Theorem 1. Note that A' 
could be not trace preserving. Rewriting Eq. ©, 

p-p(7®A s )(P+) = (l-p)(I®0)(P+) = (l-p)W. (9) 
Without loss of generality, let A s (-) = £\ . i% (•)-&],,•> 



where Ei, 



ij * 

and {l*)}^^ 1 are computational bases 



where c p is a non-negative number. ■ 

Let F = p + (1 - Cp)/. By Eq. (EL = I-F. This is 
the form of EWs in Refs. and |2|. Clearly, all EWs 
in [2H can be constructed in the form of W — p — c a I . 
A general discussion can be found in Refs. 0] and [Hj]. 
Moreover, we can also easily obtain Theorem 2 mathe- 
matically from the EW in the form P — cl (Hermitian 
matrix), where P is a positive matrix. However, they all 
have no physical interpretation. 



We can prove that A mm (p) < c p < d mm {p) if W = 
p — c p I is an EW, where A mm (p) and d min (p) are the 
minimum eigenvalue of p and the minimum diagonal ele- 
ment in p, respectively. However, for any density matrix 
p, such as the diagonal state (its density matrix is a di- 
agonal matrix), the EW in the form W = p — c p I does 
not always exist. Generally, it is not easy for any density 
matrix p to find c p to make W = p — 
have the following operational result 



c p I an EW, but we 



Theorem 3. A Hermitian matrix W = p — c p I is an 

EW for any density matrix p, if A mm (p) < c p < c™ ax = 

EyMil'l/ilV^-E^aCldilHI/.fl/d^IPe^)!- 

Ei< fcJ ;2(K| 2 |4| 2 + |/j| a |/l| 2 )|iteO*i«)|, where 

Ef 4 l^l 2 = 1, Ej B l/jf = 1 and Re(p m ) is the real 
part of the element piju of p. 



Since any density matrix p can be written 
as 9 = T, ijk i a ijki + Hyna&kb&lbijkli, where 
dijki and bijki are real, b ijk i = -b kU] , and 
tr[{\p A v B ){pAVB\)Y< l3 ki, ll tkk 31 ti b ijki\ = 0, we can 
only consider real parts of a density matrix. 



Proof. Any density matrix p on Hab can be defined 



p=52(ij\ P \M)(\m®\j)(i\) 

ijkl 



(ii) 



by computational (real orthonormal) bases {|«)}f=o 1 an( l 
{Jfc)}^- 1 in Ha, and {li}}^" 1 and {lOfe 1 m H B . 

For any unit product vector \pa^b) = Ey ^»/jK)b') 

d.A 1 j 12 1 ] v^df 



on Has with |^| 2 = 1 and Ef l/.f 



1, by Eq. 



3 



CD, 

tx\p{\HAVB){^AVB\)\ 

= (pavb\p\pavb) 

ijkl 

= E 1Re{d*f*d k fi)Re(p ijH ) 

i<k,jl 
i—k,j<l 

+ E Re(d*f;d k fi)Re( Pm ) 

i—k,j—l 



(12) 
(13) 

(14) 



(15) 



£ 2[Re(d*/; ( 4/ ; ) + M*I 2 I4| 2 + l/jfl//! 2 - Kl 2 |4| 2 



i<k,jl 



-\f J \ 2 \fi\ 2 ]Re( Pllkl )+ E 2Md*f]difi)MPw) 

i,j<i 



> E W 2 I//a™ - E 2 (l^l 4 + \fj\ 2 \fi\ 2 )\Mpm)\ 

ij i,j<l 



E WIV 

i<k 7 jl 



k\ 2 + \fj\ 2 \fi\ 2 )\Mpm)\, 



(17) 



where p^y are the diagonal elements of p. Thus, if W = 
p — CpJ, tr(VF cr) > for all separable states a.M 

A general procedure of detecting entangled states. Gen- 
erally, for any entangled state n, it is not easy to find p to 
make W — p — c p I detecting it, but we have the follow- 
ing operational result. Let us recall a well-known trace 
inequality for Hermitian matrices. 

Lemma 3. ( [26l. l27l| ) For any two Hermitian matrices 
H, K in Hab, 

d AB — l cL ab— 1 
E \(H)X dAB -i(K) < tr(HK) < E K{H)K{K), 

i=0 i=0 

where \ X {H) > ••• > X dAB . x {H), \ X {K) > ••• > 
^d AB -i(K), and Xi(H) and \i{K) are the eigenvalues 
of H and K, respectively. 

If W = p — c p I is the EW for any entangled den- 
sity matrix tt, tr(M / 7r) = tr(p7r) — c p < 0. It requires 
tr(p7r) as small as possible and c p as big as possible to 
make tr(p7r) — c p < 0. By Lemma 3, it is not diffi- 
cult to conclude that the minimum of tr(p7r) is equal 

to X)i=o _1 Xi(p)Xd AB ~i{ 7r ) if an d only if p and tt are si- 
multaneously diagonalizable [26l- 28 1 . Therefore, we have 
the following result. 

Theorem 4. An EW can be built up as commuting 
with a given entangled state. 

By Theorem 4, we have a general procedure of con- 
structing EW for any density matrix. 



(i) Suppose the spectral decomposition of any density 
matrix n = £\ Xi\ipi)(ipi\ with A < Ai < • • ■ < A^-i- 
(ii) Construct p = J2i^d AB -i-i\ipi)(^i\ with < 7 < 
7i < • • • < ld AB -i different from A < \\ < ■ • ■ < 
Xd AB -i- (hi) Compute c™^ by Theorem 3. (iv) Let 
W = p - c™ ax 7 be the EW. (v) If tr(Wn) < 0, tt is 
entangled, otherwise repeat (ii). 

Let us consider a simple example. Suppose tt p = 
pl^l + (1 -p)//4, where |V>> = ^(|00) + |11)) and 
< p < 1. It is well known that 7r p is an en- 
tangled state if p > i with positive partial transpo- 
sition (PPT) criterion^)]. Let us detect n p with the 
general procedure above, (i) The spectral decomposi- 
tion 7Tp = ^E| Wl )( Wl | + VM(w 2 | + ^|«3><w 3 | + 
VWHI. where = ^{1,0,0,-1},|W2) - 

^{1,0,0,1» 3 ) = {0,1,0,0» 4 ) = {0,0,1,0}. 
(ii) Construct p q = i+32| Wl )( Wl | + ^|wa)<a>a| + 



— 2\u 3 ) (u> 3 \ + -^{u^} (lo 4 \ and — | < g < 0. It is interest 
ing that {1^,-^, 1^,-^} >- {I 



for p > -i and 



4 ' 4 



3g 



} 



majorizes y [26|-[28j. (iii) We can compute c™' 



4 ' 4 ' 4 ' 4 

< q < 0, where "x >- y" means x 

max 1+g 

~ 4 

by Theorem 3 (see Appendix), (iv) Let W = p q — ^4^7. 

(v) tr (Wp p ) = (3p ~ 1)g < if and only if p > ± and 
— -i < g < 0. Therefore, 7r p is entangled if and only if 
p > |. Let c Pg = Interestingly, VF = p q - 

can detect w p for p > | . 

Note that we need to construct eigenvectors besides 
eigenvalues of p if there exists one or more zero eigenval- 
ues of 7r , because only if the multiplicity of each eigen- 
value of 7r is one, the spectral decomposition of 7r will 
be unique. For example, Bell state \rp) = ^75 (1 00) + 

|11)). Its density matrix Pf = \ip)(i[i\. Construct 
p = a\ip){ip\ + &[01><01j + 6|10)(10| + b\(j>){4>\, where 
\4>) = -t=(-|00) + |11)), b > a > and a + 36 = 1. 

We can compute W = p - ^1 and tr( WPf) < 0. 

Note that c p may not exist for the density matrix p 
and W = p — c p I may not be the entanglement witness 
that can detect the given n. The procedure never stops 
if the given state is not entangled (separable). By con- 
sidering a countable set of product vectors spanning the 
range of a given state, Hulpke and Brufi[29] construct an 
algorithm (procedure) to detect the separability of the 
state, in which its termination is guaranteed if the state 
is separable. Similar to [29j |. we can run both the pro- 
cedure above and the Hulpke and Brufi[29] algorithm for 
the detection of a separable state in parallel with detect- 
ing entanglement. A general discussion can be found in 
Refs. [EGA Hi- 
lt should be stressed that there is no universal EW 
that detects all entangled states, and there is no general 
procedure for constructing EWs [12f . Many problems on 
the procedure above are under investigation. 

The necessary and sufficient condition in theory. What 
is the necessary and sufficient condition for W = p — c p I 
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to be an EW in theory? Let us start from the following 
results. 

Theorem 5. For any density matrix p with spectral 
decomposition p = J2 r X r \ip r ){tp r \ if c p < c™ ax for any 
unit product vector \ha^b), tt(Wa) > for all sepa- 
rable states cr, where c™ ax = inf^^x^^i Y, r K \\ 
(VvIma^b) || 2 and W = p - c p I. 

Proof. For any unit product vector \havb), by Eq. 
and p — J2 r K\ipr)(tpr\, we have 



tr[W(\HAVB){HAVB\j\ 

(pavb\p\pavb) - c p 

(HAVb\C^2 -V#r)(^r|)|/M^s) 
r 

^] K || (VvI/M^b) IP "Cp 



> _max 
- L P 



C p > 0, 



(18) 
(19) 

(20) 
(21) 



where c" 



= "l/||p A || = l : ||,/ B || = l Er A 'r II (iPAfJ-AVB) IP- 

p — CpI, tr(W<j) > for all separable states 



c p I 



max 
P 



Thus, if W 
■ 

Corollary 1. For any density matrix p, W — p 
is an EW if and only if A min (p) < c p < 
A mm (p) is the the minimum eigenvalue and c 

W/||Mit||=l,lks||=lEr^ II (^APAVb) || 2 - 

Before giving our result in theory, we need the follow- 
ing lemmas. 

Lemma 4. ([30]) A linear map 0: Ha — > 1~Lb is posi- 
tive if and only if there exists Co,...,Cfe_i, I?o,---,-D/-i on 
T~Lab such that {Dj}jl} are a contractive locally linear 
fc-i 



combination of {Ci}f =0 and 



fe-i 



(22) 



i=0 j=0 

for all X on T-Lb- 

Lemma 5. [30] Furthermore, O in Eq. 
pletely positive if and only if {Dj} l ~} is a linear combi- 
nation of {C'i}i~Q with a contractive coefficient matrix. 

Theorem 6. A Hermitian matrix W — p — c p I is 
an EW for any density matrix p = ( I ® A) (P+ ) if and 
only if {-y/ d-ABCpE, 



|* }{_ o is a contractive locally linear 



combination of {J7 r }^ZQ but not a linear combination 



of {J7 r } r _o with a contractive coefficient matrix, where 



A(- 



and 



computational bases in Ha, {\j}}i=o 1 m ^s- 
Proof. Any density matrix can be written as 

p = ^2 A r |Vv)(Vv| 



(23) 



by means of its spectral decomposition with nonnegative 
eigenvalues A r . Taking V r such that J® V r \/3) = \tp r ), 
any density matrix p can be written as 



where A(-) = J2 r U r (-)U~£ is a completely positive map 
and U r = \fX~rV r , I®V r \f3) = \%p r ), (n\V r \m) = y/dABOmn 
and |-0 r ) = J2 mn Omn|m)|n)5i|. 

By Eq. ©, VF = I <g> A(P+) - J ® A' a (P + ), where 
Ai(-) = Et v/^^^OXyS^^JV- % Lemmas 
4 and 5, { \J (Iab CpE^f }jZ is a contractive locally linear 
combination of {U r }rZo but not a linear combination of 
{U r } with a contractive coefficient matrix if and only if 
$ = A - A' s is a PNCP map, and W = p - c p J is an EW. 
■ 

In conclusion, we have demonstrated that any EW can 
be constructed from a certain density matrix. This re- 
sult shows that the research on density matrices can re- 
place the research on entanglement witnesses. The trace 
inequality reveals the general procedure of constructing 
EW for any density matrix. Both in operation and in 
theory, the necessary and sufficient condition of an EW 
in the form W = p — c p I and some examples are given. 
Here we only consider the bipartite case on the finite di- 
mensional Hilbert space, but we can also generalize our 
results to multipartite system and infinite dimensional 
Hilbert space. 
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Appendix: The Procedure of Computing c 



I'; 



Clearly, 



tr(p q (\fi A )v B ) (li a (v b \)) = (a* A {^B | Pq | Pa) v B ) 

= k/o(ooi +d$/r<oii +d 1 / *(ioi +d;/r<ni) x Pq x 

(do/o|00> + do/i|01> + rfi/o|10) + di/i|ll» 



+Mi| 2 l/i 



4 

,1 + q 



IrfilU 



i + d* fo di fx | + d\ fi d f - 



(26) 



(27) 



|do| 2 (|/o| 2 + |/i| 2 ) 



l + q 



+ Mi| 2 (l/o| 2 + |/i| 2 ) 



4 

l+q 



J0 | 2 |/i| 2 2q 

Mi| 2 l/o 



4 



2Re(dSdi/ */i) 



^(Mop + Mrnd/oP + IAI 2 )^ 

= i+ff + [R e(d-/ 1 )-He(d 1 /o*)] 2 ^ 
1 + d 



[Re(d* h) - Re(difo)] 



> 
~ 4 

Therefore, we have c™ ax = — r 2 . In addition, we have 
another method for computing c™ ax . 

Any qubit pure state \ip) can be written as \ip) = a|0) + 
(3\1), where a and (3 are complex number and |cy| 2 + \f3\ 2 = 
1. Because |a| 2 + |/3| 2 = 1, \yj) can be rewritten as 



2 

2-JL 
2 

(28) 
(29) 



|V,)=e"\cos-|0> + e-sin-|l)), 



(30) 



where 0, r and s are real numbers. The factor of e lr out 
the front can be ignored since it has no observable effects 
[iH, and for that reason, \ip) can be effectively written 
as 

| V , )=cos ^|0)+ e - sh/|l). (31) 
Therefore, \pa)vb) for two qubits can be written as 



|MAW> = (cos^|0) + e iSl sinJ-|l))(cos^|0)+e^ sin^|l)) 



Pq 



/ 1+2 

/ 4 




£ 

2 





V 

o V 





2 

2 




1+2 
4 



where — | < q < 0. 



(25) 



By Theorem 3, \pa)vb) for two qubits can be written 
as \va)vb) = do/o|00) + d /i|01) + di/ |10) + di/i|H) 
with |do| 2 + |di| 2 = 1, |/o| 2 + |/i| 2 = 1. 



01 °2 lnn\ 

cos — cos — 1 00) 

. 0\ 6 
+e 1 sin — cos - 



e lS2 cos -i sin |01> 



2 2 
2 |10) +e 4(si+S2) 



(32) 



sin y sin y |11) 



to[p q (\p A )vB){p-A{vB\). 
l+q 2 6 x 2 9 2 
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